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There will be two interesting panel discussions at TIME 2012: 

DERIVE is Still Special 

Organizers: Michel Beaudin, Josef Bohm, Jose Luis Galan, Carl Leinbach 

Abstract 

At the TIME 2006 conference in Dresden, several DERIVE program users were shocked to hear 
that DERIVE would be discontinued and replaced by the new product TI-Nspire CAS. In 2007, a 
transition period started where users asked TI developers to implement some interesting features 
of DERIVE. TI-Nspire CAS, with two platforms and a wider spectrum - due mainly to numer- 
ous integrated applications-, could not satisfy all wishes but tried to address some of them: one 
particular example is found in the latest operating system of Nspire CAS (OS 3.2) where 3D pa- 
rametric plotting has been added with David Parker on board. But there are still special and 
unique features of DERIVE that cannot be found in Nspire CAS. 

In this panel discussion, we would like participants attending this session to talk about what they 
still find special in DERIVE: to do so, live examples in various fields of mathematics will be 
shown to the audience using DERIVE. One goal of this discussion is to provide constructive 
ideas to Texas Instruments developers: proclaimed at its beginning as the “true successor of DE- 
RIVE”, Nspire CAS has all the potential needed to satisfy many requests of users. 

Topics for discussion include, but are not limited to: Subexpressions, Showing Step Simplifica- 
tion, 2D implicit plotting, precision digits, cubic equations, piecewise continuous functions, Ex- 
pand/Collect for trig expressions, creating random numbers without the randseed, plotting poly- 
hedrons in 3D, iteration, Fourier series, polynomial systems, Grobner bases. Those who would 
like to present one’s favourite DERIVE special which has not been implemented in TI-Nspire 
so far are welcome to join. 

I would like to ask you about your DERIVE SPECIAL which you are still missing in 
TI-NspireCAS. We know that the Nspire-developers are eager to implement as many DERIVE 
features as possible - but this needs time. Much has been done since the first presentation of 
TI-Nspire in Dresden 2006 but the DERIVIANs are still waiting for some improvements. 

We will publish your wishes together with the TIME 2012 results in our newsletter. 


This is a second panel discussion (part of the abstract): 

Assessment and mathematical placement exams with technology 

An effective yield of the use of technology in mathematics education can only be ex- 
pected if the tools are also applied accordingly in the exam situation. Discussing this 
topic we have to keep in mind on the one hand, the process oriented assessment as an 
integral part of the learning process, and on the other hand, product oriented mathemati- 
cal placement exams and final exams as a proof of authority for the students. 

Your contributions would be very welcome , Josef. 
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Liebe DUG-Mitglieder, 

nun ist es mir doch tatsachlich gelun- 
gen, den DNL#86 noch vor der TIME 2012 
fertig zu stellen. Die revidierte Fassung 
des DNL#29 muss da fur noch etwas war- 
ten. 

Im Rahmen der TIME 2012 wird es 
zwei interessante Podiumsdiskussionen ge- 
ben: Wir werden uber DERIVE-features 
sprechen, die wir in der aktuellen TI- 
NspireCAS-Version noch vermissen. Die 
zweite Diskussion betrifft Probleme und 
Chancen von Technologie unterstiitzten 
Priifungen. Hinweise auf die Diskussionsin- 
halte finden Sie auf der Informationsseite. 
Sie sind herzlich eingeladen, Ihre Meinung 
oder Beitrage zu beiden Themen auch im 
Rahmen des DNL abzugeben. Im Herbst 
werde ich uber die Ergebnisse der Diskus- 
sionen berichten. 

Nun einige Hinweise zum Inhalt dieses 
DNL: Guido Herweyers' schones Statistik- 
skript wird weiter gefiihrt. Michel und Gil- 
les haben in Malaga einen spannenden Vor- 
trag zur numerischen Losung von Differen- 
tialgleichungen 2. Ordnung gehalten. Da sie 
sich auf die Behandlung mit dem V200 be- 
schrankten, habe ich versucht, ihre Ideen 
mit DERIVE bzw. mit TI-NspireCAS auf- 
zugreifen. 

Eine wiederum sehr gelungene schiiler- 
gerechte Untersuchung wird von Roland 
Schroder beigesteuert, wahrend der Arti- 
kel unseres kolumbianischen Freunds, Nel- 
son Urrego, wohl schon etwas weiter von 
der Schulmathematik entfernt ist. Zu sei- 
nem Thema habe ich noch ein paar sehr 
brauchbare websites gef unden. 

Dietmar Oertel hat mir eine ganze Kol- 
lektion von DERIVE-Grafiken geschickt, die 
er im Zusammenhang mit seinen Untersu- 
chungen mit Taylorreihen in stundenlangen 
Sitzungen erzeugt hat. Einen kleinen Blick 
in seine Galerie konnen Sie auf Seite 42 
werfen. 

Ich wiinsche einen schonen Sommer und 
freue mich auf den DNL#87. 


Dear DUG Members, 

yes, I was successful in finishing this 
DNL before start of TIME 2012. The re- 
vised version on DNL#29 must wait a little 
while. 

In the frame of TIME 2012 we will 
have two interesting panel discussions: we 
will speak about DERIVE-features which we 
have been missing in the NspireCAS ver- 
sions until now. The second discussion will 
focus on the problems and chances of tech- 
nology supported assessments. Abstracts 
of both discussions can be found on the 
information page. You are invited to con- 
tribute or to comment both subjects in our 
DNL. I will report about the results of the 
discussions an publish your possible reac- 
tions in the fall-issue of the DNL. 

I have some remarks on the contents of 
this Newsletter: You can find the second 
part of Guido Herweyers’ excellent statis- 
tics paper. Michel and Gilles gave an excit- 
ing lecture on the numerical solution of 2 nd 
order ODEs in Malaga. As they showed the 
V200-treatment only I tried to follow their 
ideas with DERIVE and TI-NspireCAS. 

Roland Schroder contributed again a 
very fine student centred investigation. 
Our friend from Colombia, Nelson Urrego 
sent a paper which has some distance from 
school mathematics. I found a couple of 
interesting websites referring to the con- 
tents of Nelson's article. 

Dietmar Oertel sent a collection of 
DERIVE graphs produced by him in connec- 
tion with his investigations of Taylor series. 
I believe his words when writing that he 
needed hours and hours plotting the fami- 
lies of curves. You may take a tiny look on 
his gallery on page 42. 

Finally I’d like to wish you a fine sum- 
mer and am looking forward to presenting 
DNL#87 in three months 

Viele GruBe, kindest regards 


Download all DTVL-Derive- and TI-f!les from 

http : / / www . austromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 7/ User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla, Austria 
Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


I II 

Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS- 77 Newsletter will be. 

Next issue: September 2012 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Liike-Rosendahl, GER 
Hill-Encription, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the TI, M. R. Phillips, USA 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Structures in the Set of Prime Numbers, D. Oertel, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Laplace Transforms, ODEs and CAS, G. Piccard & Ch. Trottier, CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Cesar Multiplication, G. Schodl, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 
Rational Hooks, J. Lechner, AUT 
Mathematical Model for Snail Shells (4), P. Trebisz, GER 
Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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Statistics with Tl-Nspire 3.1 (Part 2) 

Visualising and Simulating Dynamically with Tl-Nspire 3.1 

Guido Herweyers, KHBO Campus Oostende 
guido.herwevers@khbo.be 


(3) Quantitative ungrouped data 

Example 4: 

Given are the following birth weights of 16 girls and 14 boys (in kg). Data from an Excel data stock 
can be imported directly into a Tl-Nspire spreadsheet page by Copy and Paste. For some countries is 
one obstacle to overcome: the decimal comma must be changed to a decimal point first. 

(From German: Extras > Optionen > International > . . .) 


3 

Boys 

Girls 

Boys 

3,54 

3,58 

3.54 

3.58 

3,49 

3,59 

3.49 

3.59 

2,72 

3,6 

2.72 

3.6 

4,13 

3,61 

4.13 

3.61 

3,58 

3,62 

3.58 

3.62 

3,36 

3,63 

3.36 

3.63 

3,67 

3,64 

3.67 

3.64 

3,22 

3,65 

3.22 

3.65 

3,22 

3,66 

3.22 

3.66 

3,22 

3,67 

3.22 

3.67 

4,08 

3,68 

4.08 

3.68 

4,58 

3,69 

4.58 

3.69 

3,13 

3,7 

3.13 

3.7 

3,44 

3,71 

3.44 

3.71 

3,4 


3.4 


1,63 


1.63 



—> Start a new problem within the same file 
(in the document menu Insert Problem), 
change the problem name Problem 2 to 

Quantitative Variables, and insert a Lists 
& Spreadsheet page. 


Copy the data into columns A and B. (Take 
care to highlight the according region of 
cells =16x2 before pasting.) 

Name the columns as girls and boys. 


Content 


Documents 






Split the screen in three windows as shown 
below and insert two Data & Statistics 
applications into the empty windows. 
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Insert - 

© - 



* girls 

B boys 

♦ 

1 

i 

1 

3.54 

3.5 

2 

3.49 

3.5! 

3 

2.72 

3d 

4 

4.13 

3.6 


m 


m 


v girls B boys c 


1 3.54 3.58 

2 3.49 3.59 

3 2.72 3.6 

4 4.13 3.61 

5 3.58 3.62 

6 3.36 3.63 

7 3.67 3.64 

8 3.22 3.65 

9 3.22 3.66 


10 

3.22 

3.67 

11 

4.08 

3.68 

12 

4.58 

3.69 

13 

3.13 

3.7 

14 

3.44 

3.71 

L® 

J boys 


Click here to add an application 


1 :Add Calculator 
£ 2:Add Graphs 
l3 3:Add Geometry 
J 4:Add Lists & Spreadsheet 


ft 5:Add Data & Statistics 


3 6:Add Notes 

7:Add Vernier DataQuest 


Click here to add an application 


—> Click in both windows on the bottom for adding the variable girls. By default appears a Dot Plot in 
both windows. Change the diagram in the upper window in a Box Plot (right mouse click = rmcl). 



— » Enter “Girls” in cell C3 of the spreadsheet (under quotes), then “Mean” in cell C4 and “Median” in 
cell C5. 

Type =mean(girls) in cell D4 followed by pressing the ( ENTER) -key; the mean appears. 

Type =median(girls) in cell D5 followed by pressing the ( ENTER -key; the median appears. 

You can do the same for the boys’ data in cells C7 to D9. 
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A girls 

B boys 

c 

D 

♦ 





1 

3.54 

3.58 



2 

3.49 

3.59 



3 

4 

2.72 

4.13 

3.6 

3.61 

Girls 

Mean 

3.40063 

5 

3.58 

3.62 

Median 

3.42 

6 

7 

3.36 

3.67 

3.63 

3.64 

Boys 


8 

3.22 

3.65 

Mean 

3.645 

9 

3.22 

3.66 

Median 

3.645 

10 

3.22 

3.67 




vl := mean(girls) 
= 3.40062 


1.4 1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6 
girts 


Documents Toolbox 

| 1 0 ■> ri 


Data & Statistics 


|5f j: Plot Type 
: 2: Plot Properties 


A girls B boys 


4: Analyze ► 

1 : Remove 

5:Window/Zoom ► 

2:Add Movable Line 


3:Lock Intercept at Zero 


3.58 


3.59 
3*6 


—> Click into the Box Plot window 
and then choose first the Analyze 
menu in the Documents Toolbox 
and then the option Plot Value. 

Fill in for vl:= mean(girls). A ver- 
tical line in the position of the mean 
will appear (right above). 

Now you can grip one outlier (left 
mouse click on it) and move it 
along. You can observe how the 
data changes together with a change 
of the mean. Investigate the influ- 
ence of the outlier on mean and 
median. 

The data are altered accordingly in 
the spreadsheet and dot plot, too. 


4: Plot Function 


5:Shade Under Function 

3.61 

' n '- u - 6: Regression ► 

3.62 

V[ 7: Residuals ► 

3.63 

S:Plot Value 

3.64 

9:Show Normal PDF 

3.65 

A:Graph Trace 

3.66 


1 A girls 

B boys 

■ 

■ ip 

1 


♦ 





vl := meanlgirlsj 


094 Q 

T 

4 

2.72 

4.13 

3.6 

3.61 

Girls 

Mean 

3.49213 

= 3.49213 

& 

5 

3.58 

3.62 

Median 

3.42 




6 

7 

3.36 

3.67 

3.63 

3.64 

Boys 



■ 

t 

8 

3.22 

3.65 

Mean 

3.645 

1! 


n — • — i — ■ — i — ■ — i — • — i — ■ — i — ■ — i — ■ — i — ■ — i — 

1.4 1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6 
girls 

9 

3.22 

3.66 

Median 

3.645 

10 

3.22 

3.67 




11 

12 

4.08 

4.58 

3.68 

3.69 




13 

3.13 

3.7 




14 

15 

3.44 

3.4 

3.71 



°ooo 

O CD 

O OOOOO G> O 

16 

3.094 




V 




5 

' ' ' '''''' 

1 * 1 ' 1 * 1 ^ 

3.4 3.8 4.2 4.6 

5 

A1 =3.54 



4 ► 

1 1.4 1.8 2.2 2.6 3.0 
girl: 
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Hint: It is easy to come back to the original data by pressing the Undo-key: 



Fortunately it is possible to forbid (to lock) altering data. Open a Calculator page and type 
lock(girls). Now try again moving the outlier. Command unlock(girls) gives the data free again. 


— » Restore the original data. Click on the 
vertical mean line and press ( DbI ) . Make 
a right mouse click on the girls below the 
horizontal axis to achieve Add X vari- 
able and then click on the offered variable 
boys. Do the same in the Box Plot win- 
dow. Adapt the colours. Compare the data 
and their representations of girls and boys 
for these samples. 



—> Now we will find all statistics of the given data. Click in one cell of column D in the spreadsheet, 
choose the Statistics menu, continue with Stat Calculations and choose option One- Variable Statis- 
tics. 
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1 1 


F 

1 J 

G 

■ 


♦ 

i 

|=OneVar('girls.1)| 

=OneVar('boys.1) 

00 

i • ■ • 

1 

Title 

One-Variable S... 

One-Variable S... 


2 

X 

3.40063 

3.645 


3 

lx 

54.41 

51.03 

| 1 

4 

lx 2 

191.317 

186.027 

-§ I 



0.647523 




5 

SX := Sn-.„ 

U. U 4 1 ojj 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 

1.4 1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6 

6 

OX := On... 

0.626962 

0.040311 

7 

n 

16. 

14. 

•girls f) boys 

8 

MinX 

1.63 

3.58 


9 

QiX 

3.22 

3.61 

girls 

10 

MedianX... 

3.42 

3.645 

• JJfi t • 

11 

Q 3 X 

3.625 

3.68 

12 

MaxX 

4.58 

3.71 

| 1 

13 

SSX := I... 

6.28929 

0.02275 


14 

15 









g 

i ■ i ■ i ■ i ■ i • i • i ■ i • i • i| 

1.4 1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6 

G 

=OneVarl, 'girls, l): CopyVar Stat., Statl. 


•girls# boys 


— » Change again the data (by dragging) and observe the changes of the statistics. 


—> Open a Calculator App and type statl. and select one statistic of 
the appearing list, followed by ( ENTER) . Do the same for stat2. in order 
to get the statistics of the second list (the boys). 

statl. results delivers a matrix giving all statistics for the girls’ data. 
TI-Nspire stores automatically all statistics calculations in statl, stat2, 
stat3, ... 


Lock girls 


Unlock girls 
stall .MedianX 


stat2 


% maxx 


medianx 


[^1 minx 
|%| n 
% qix 

q 3X 




0 
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& 

statl .MedianX 3.42 


stat2.MedianX 


3.645 

statL results 

"Title" 

"One- Variable Statistics" 



"x" 

3.40063 



"Ex" 

54.41 



"Ex 2 " 

191.317 



"sx := sn-ix" 

0.647523 



"ax := anx" 

0.626962 



"n" 

16. 



"MinX" 

1.63 



"QiX" 

3.22 



"MedianX" 

3.42 



"QsX" 

3.625 



" MaxX " 

4.58 



"SSX := E(x-x) 2 " 

6.28929 



Statistics calculations can also be performed in a Calculator App via the Statistics menu. 


Example 5: 

In this example ungrouped quantitative data are split according categories. The population is the data 
set of the first year bachelor students at the KHBO Campus in Oostende in the academic year 2008/09 
(137 students). It follows an overview of the results of a mathematics test in the first semester (maxi- 
mum score 10 points) together with the type of secondary school attended by the students. There are 
two lists: results and school_type. 


|P results 

® school type 

* 1 

i i 

i' 

3 

Chemie TSO 

2 

5 

Wet-wisk. ASO 

3 

5.5 

Wet-wisk. ASO 

4 

5 

Wet-wisk. ASO 

5 

6 

Latijn-wisk. ASO 

6 

7 

Latijn-wet. ASO 

7 

2Techn. wet. TSO 

8 

2 

Chemie TSO 

9 

7 

Wet-sport ASO 

1 

4 

Wet-wisk. ASO 

ii 

4.5 

Humane wet. ASO 

12 

4.5Techn. wet. TSO 

13 

7Talen-wet. ASO 

1 

4.5 

Wet-wisk. ASO 

15 

4 

Biotechn. wet. TSO 

B 

schooljype 


Insert a new Data & Statistics page. 


D 
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Add variable results for the x-axis and variable school_type for the y- axis. This gives a Dot Plot of 
the test results split according to the various school types. Now add a vertical line for score 5 (again in 
Documents Toolbox > Analyze > Plot value and change the variable name vl to score). 


Biotechn. wet. TSO 
Buitenland 



O | O O r 

I score:=5 

O I o 

Chemie TSO 


O 

o 0 1 

Econ-talen ASO 



I o 

El.inst.techn.TSO 

C 

) 0 

0 0 O O I 

Elektr.mech. TSO 


O 

0 O O 0 4* 0 

Elektri-elektro TSO 



o 

c 

o 

o 

o 

Humane wet. ASO 



O I 

J IW TSO 



| o o o o o o 

§ Koel-warm. TSO 



° ! 

-g Latijn-wet. ASO 



o 1 o 

Latijn-wisk. ASO 



o I o o o o 

Mech.vormg. TSO 



O O I 

Talen-wet. ASO 



o o o ] o o 

Techn. wet. TSO 


o 

o 0 o $ ooooo 

Vliegtuigtechn. TSO 


o 


Wet-sport ASO 



I o 

Wet-wisk. ASO 



o o o 4> o o OOO oo 


c 

1 

) 1 

— 1 1 1 T — 1 1 1 1 1 

23456789 10 




results 


The dot plot gives a clear impression of the individual results and which secondary school types 
reached low or high scores. 

—> Change the Dot Plot into a Box Plot. We can read off that there are six types with a median greater 
than score 5 and 12 types with a median less 5. 


Biotechn. wet. TSO 

J score:=5 

Buitenland 


Chemie TSO 

— I 

Econ-talen ASO 

1 < 

El.inst.techn.TSO 

-i ■ i I 

Elektr.mech. TSO 

1 1 A — 

g_ Elektri-elektro TSO 

— 1 — ■ — j 

& Humane wet. ASO 

> i 

-i IW TSO 

1 — 

§ Koel-warm. TSO 

1 | 

-g Latijn-wet. ASO 


vi Latijn-wisk. ASO 

1 1 1 

Mech.vormg. TSO 

CD | 

Talen-wet. ASO 


Techn. wet. TSO 

U 1 j 1 

Vliegtuigtechn. TSO 

1 I 

Wet-sport ASO 
W6t — wisk ASO 

1 1 

C 

1 1 1 1 1 1 

)1 23456789 10 


results 


p 10 
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Remove the variable results from the x-axis (rmcl on results). 


Elektri-elektro TS 
Elektr.mech. TS 
Techn. wet. TS 
Wet-wisk. AS 
El.inst. techn. TS 
IWTSOI 
Talen-wet. AS 

<L> 

Q- Latijn-wisk. ASC 
.JBiotechn. wet. TS( 
? Chemie TS( 


Buitenlan 
Latijn-wet. ASO l | 2 (1 .5%) 

Mech.vormg. TS d | 2 (1 .5%) 

Econ-talen AScQl (0.7%) 
Humane wet. AScj |l (0.7%) 
Koel-warm. TSd |l (0.7%) 
Vliegtuigtechn. TSd |l (0.7%) 
Wet-sport AScBl (0.7%) 



5 (3.6%) 
4 (2.9%) 

4 (2.9%) 

2 (1 .5%) 


T 

2 


~r~ 

4 


— i — 

10 


— i 

1 2 


— i 1 1 — 

14 16 18 

Frequency 


— i — 
20 


— i — 
22 


— i — 

24 


— i — 
26 


— i — 

28 


— » Insert a new page and split it horizontally into a Data & Statistics App and a Calculator App. 
Collect all data in one single Box Plot and in the statistics of the 137 data records. 



OneVar results, 1 : stat. results 



"sx := Sn-ix" 

2.21038 

"ax := anx" 

2.2023 

"n" 

137. 

"MinX" 

0 . 

" QiX" 

3. 

"MedianX" 

4.5 

"Q 3 X" 

6 . 

"MaxX" 

9.5 

'SSX := Z(x-x) 2 " 

664.467 


1/99 


i 1 1 1 1 1 1 1 1 1 r~ 

01 23456789 10 

results 


’Title" 

"One-Variable Statistics 

"x" 

4.66058 

"Zx" 

638.5 

"Zx 2 " 

3640.25 
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Example 6: 

During summer holidays the Belgian Coast is watched by life guards. Many students feel attracted by 
this adventurous holiday job. Apprenticeship for life guard is not to be underestimated. In addition to 
an extended theoretical course a couple of difficult swimming tests must be passed. 

The training for life guards is organized annually. At first the students must pass the theoretical exam. 
The course exists of seven main parts which are all assessed separately - each of them with 100 points 
maximum. For passing the theoretical part at least 50% of each main part must be reached. The stu- 
dents who have passed the first part can go on with the practical part of the exam (swimming tests, 
knowledge about knots, first aid in case of accidents). The results of the theoretical exam are published 
on a website annually [4]. 

There were 367 participants for the course in 2010-2011. Each of them got a number between 1 and 
367. The results are given sorted wrt the numbers in the TI-N spire-file. 

— » Insert a Calculator page and lock the data with the command 

lock hl,h2,h3,h4,h5,h6,h7 

—> Insert a new page with two Data & Statistics applications. Produce a Dot Plot and a Box Plot as 
well of the results of main part 4 (given in variable h4). 



—> Change the Box Plot to a Histogram, enlarge the histogram window (by dragging the horizontal 
borderline). Experiment with the class width (by dragging a vertical side of any rectangle). One can set 
the class width manually, too (rmcl on the histogram, Option 5: Bin Settings). Move the cursor over 
the histogram in order to read off the frequencies. You can represent the relative frequencies making 
again a right mouse click on the histogram and choosing Scale > Percent. 

When selecting one class the corresponding data are marked in the dot plot. Clicking in the free space 
outside of the diagram cancels the selection. 
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—> The class frequencies can be calculated in a Calculator page using the command countlf. You will 
obtain the frequency list for the classes [0,10[, [10,20[, . . [100,1 10[ with the command line 
seq(countIf(h4, 10k < ? < 20(k+l), k, 0, 10). The list can also be generated in a column of the spread- 
sheet by typing =seq(countIf(h4, 10k < ? < 20(k+l), k, 0, 10) into the grey cell in the second row of 
the column. 


Lock hi, Ji2, 1i3,1i4Ji5,1i6 f Jt7 

w 

Done n 

countlf(//4,20<?<30) 

11 

seq(countIf(A4,10-ifc<?<10-(fc+l)),lt,0,10) 

{ 2,7, 1 1, 14,27,58,70, 100,67, 1 1,0 } 

1 


— » Compare the results for all main parts in one window by their box plots and discuss the differences. 



ro 

_c 



_c 


in 

-C 


ID 

-C 


m m •• 




lx 

-C 


0 



i 1 1 1 1 1 1 r - 

10 20 30 40 50 60 70 80 

»h1 #h2Qh3#h4#h5Qh6 #h7 


“i 1 — 

90 100 
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—> Call the statistics of the results of first three main parts of the exam. Find the connections to the 
respective box plots. 


OneVar 3,/z7,/z2,/zi: stat.results 





"Title" 

"One-Variable Statistics" 

11 11 

H 11 



"X" 

59.9019 

81.0198 

57.8392 



"Ex" 

21984. 

28681. 

21227. 



"lx 2 " 

1.43359e6 

2.36305e6 

1.38082e6 



"sx := Sn-ix" 

17.8566 

10.5536 

20.4505 



"ox := anx" 

17.8323 

10.5387 

20.4226 



"n" 

367. 

354. 

367. 



"MinX" 

12. 

33. 

0 . 



"QiX" 

50. 

74. 

43. 



"MedianX" 

60. 

82. 

62. 



"Q3X" 

73. 

89. 

74. 



"MaxX" 

100. 

100. 

99. 



"SSX := Z(x-x) 2 " 

116702. 

39316.9 

153070. . 



Lists h2 and h7 contain empty cells for the students who didn’t participate in the exams of these main 
parts. 

Empty cells of a list (or of a column in a speadsheet table) are marked by an underscore ( _ ) or by the 
word void. 


dim|fi2) 

A 

367 ' 

meanU/2) 

28681 

354 

sum(ft2j 

28681 

367 

367 

isVoid(ft2) 

{ false, false, false, false, false, false, false, false, false, false, false, false, false, false* m 

countIf|isVoid\ft2J,?=false) 

354 

1 ' 1 'i 'i 

c o u ntl f| i s V o i d\ ft2 J , ? =tru e J 

13 

li2 

{%, 97, 78, 83, 90, 85, 85, 92, 68, 90, 86, 63, 93, 78, 97, 95, 88, 79, 89, 93, 95, 79, 66 , 89 


There are 13 empty cells in list h2. 
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(4) Quantitative grouped data 

Example 7: 

Given are the frequencies for main part 4 results for the classes as shown in the table as follows: 


[0,10[ 

[10,20[ 

[20,30[ 

[30,40[ 

[40,50[ 

[50,60[ 

[60,70[ 

[70,80[ 

[80,90[ 

[90,100[ 

2 

7 

11 

14 

27 

58 

70 

100 

67 

11 


—> Create a histogram based on this frequency table. Find the relative frequencies and the cumulative 
frequencies as well. Draw the relative frequency graph. 



® Left cutpts 

B midpts ®freq 

D ret freq 

E cum ret freq 

* 

®seq(1O*k,kA10) 

= seq ( co u ntif (' h 4 , 1 0*k 


=round £freq/(sum (freq)),3) 


1 

0 

5 

2 

0.005 

0 

2 

10 

15 

7 

0.019 

0.005 

3 

20 

25 

11 


0.03 

0.024 

4 

30 

35 

14 


0.038 

0.054 

5 

40 

45 

27 


0.074 

0.092 

6 

50 

55 

58 


0.158 

0.166 

7 

60 

65 

70 


0.191 

0.324 

8 

70 

75 

100 

0.272 

0.515 

9 

80 

85 

67 

0.183 

0.787 

10 

90 

95 

11 

0.03 

0.97 

11 

100 



1. 


The second column contains the list class of the midpoints of the class intervals. This list is easy 
done: click in cell B1 and type 

= (al + a2)/2 

followed by ( ENTER) (start a formula with an equation sign like in MS-Excel). Now you can copy 
down this formula until cell BIO. For doing this move the cursor to the right bottom corner of cell B1 
until a plus-sign is appearing. Keep the left mouse key pressed and move the cursor down to BIO to 
spread the formula. 

Column E is filled with the cumulate frequencies: write 0 into cell El and the formula = el+dl in cell 
E2. Copy down the cell contents until cell E10. 
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(5) Scatter diagrams and regression 

Example 8: 

Given is the table containing the percentage of Belgian homes having internet access (source NIS): 


Year 

2005 

2006 

2007 

2008 

2009 

Percentage 

50 

54 

60 

64 

67 


—> Represent the given data in a scatter diagram. 

Open a new Lists & Spreadsheet page and fill in the data (year = 0 corresponds with year = 2005). 



A year 

B percent... 

♦ 

1 



0 

50 

2 

1 

54 

3 

2 

60 

4 

3 

64 

5 

4 

67 


Insert a Data & Statistics page. Set year on the horizontal axis and percentage on the vertical one. 
Add a Moveable Line via the Analyze menu. 


Data & Statistics 


[£f 1: Plot Type 
[j}| 2: Plot Properties 
1^- 3:Actions 


u> 


5:Window/Zoom 


3: Lock Intercept at Zero 

IV 4: Plot Function 

yV 5: Shade Under Function 

^^6: Regression ' 

7: Residuals ' 

8:Plot Value 

yjV 9:Show Normal PDF 

K A: Graph Trace 

54 - 


50 - 




0.0 


0.5 


1.0 


1.5 


2.0 

year 


2.5 


3.0 


3.5 


4.0 


Continue in the Analyze menu with option Residuals > 1: Show Residual Squares. 
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Move the line (You can move the line by gripping it around the midpoint of the visible segment - the 
cursor changes into a fourfold arrow - and you can rotate the line by gripping at the ends of the seg- 
ment - the cursor changes in a circle formed arrow). Play around and try to find the minimum value of 
the Sum of Squares by experimenting with the position of the line. 

Finally compare your “best straight line” and its corresponding sum of squares with TI-Nspire’s re- 
sults. Use the menu Analyze > Regression > Show Linear (mx + b) followed by Analyze > Residu- 
als > Show Residual Squares. 



Return to the Data & Statistics page (by pressing Ctrl + <— or via the page sorter). Click into one cell 
of column C and choose menu Statistics > Stat Calculations > Linear Regression (mx + b) ... 
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A year 

B percent... 

C D 

E 

■ 

♦ 

1 






0 

50 




2 

1 

54 




3 

2 

60 







4 

3 

64 





5 

4 

67 


X List: ('year 
Y List: ('percentage 
Save RegEqn to: |f1 


6 




; " | 


7 




Frequency List: |l 



8 




Category List: 1 


9 


1 

L 

Include Categories: 


10 




1st Result Column: ld[] 

”1 


11 




OK | Cancel | 





The coefficient of determination r 2 = 98,8% is excellent. The equation of the regression line is stored 
as function fl. Estimate the percentage for 2010 in cell F9 applying the formula =fl(5). 



Example 9: 

The graph of a time series is a scatter diagram with data of a quantitative variable on the vertical axis 
as function of time on the horizontal axis. The points are connected by segments. In this way the trend 
of development of the variable can be studied. 

The tables given below show the development of the numbers of weddings and divorces in Belgium in 
the years 1995 to 2009. 


1 

Aantal huwelijken per gerj 

peente, 1995 - 2009 


Reims ADM W ISTRATE VE EENHED6N 

1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 


TDtaal 

51.402 50.552 47.759 44.393 44.171 45.123 42.110 40.434 41.777 43.296 43.141 44.813 45.561 45.613 43.303 


Bron (verplichte vermeldmg ): Algemenej 

iDirectie Statistiek en Economische In formatie -Thematische Directie Samenleving. 



Aantal echtscheidingen pier gemeente, 1995 - 2009 


Refnis|AOMMSTRATEVE EENHEDEN 1995| 1996| 1997| 1998| 1999| 2000 1 2001 1 2002 1 2003 1 2004 1 2005 1 2006 1 2007 1 2008 1 2009 


-otaal 1 34,983 28.402 26.748 26.503 26.423 27.002 29.314 30.628 31.355 31.405 30.840 29.189 30.081 35.366 32.606| 

Bron (verplichte vermeldmg): AlgemenejOirectie Statistiek en Economische Informatte -Thematische Directie S amen levin g. 
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— > Represent both time series and additionally represent the differences of the data. 


A time 

B weddings 

C divorces 

D difference 

■ ■ 

■ 

■ 

a 

♦ 




=, weddings-'divorces 




1 

1 

1995 

51402 

34983 

16419 





2 

1996 

50552 

28402 

22150 






3 

1997 

47759 

26748 

21011 






4 

1998 

44393 

26503 

17890 





5 

1999 

44171 

26423 

17748 



1 



A2A15 

=al + 1 




50000 - 


140000- 


>30000 

ID 

□ 

J20000- 

I 

£10000 



• weddings 
□ divorces 
A difference 


1994 


1996 


1998 


2000 


2002 

time 


2004 


2006 


2008 


201 C 


After typing the name difference for column D automatically appears differences when clicking in 
the grey formula cell below the caption of the column. 

Accomplish the formula with 


differencesweddings - divorces 


Hints: 

The names of the variables are offered in a list via the ED key 
in document menu. Then it is easy to insert a name into a 
formula. 


■■time 

■ weddings 

■ divorces ■ difference B F 

* 


| := 'weddings -'divorces 

H 1:class 
iH;2:classborder 
H 3:cum_rel_freq 
H 4:cumrel_freq 
H 5: difference 

1 1995 

51402 

34983 16419 

2 1996 

3 1997 

50552 

47759 

28402 22150 

26748 21011 

4 1998 

44393 

26503 17890 

[ffi 6: divorces 

5 1999 

44171 

26423 17748 

fci7:echtscheidingen 
tfwl 8:f1 
H 9:freq 




D difference : = 'weddings-' divorces 


The legend does not appear 
automatically. 


1:Zoom 

2:Hide Connecting Lines 


3:Show Legend 


4: Remove weddings 
5: Remove divorces 
6: Remove difference 

^ 


You need Add Y Variable to present more than only one time 
series. 


Next DNL: Part 2: Simulations with TI-NspireCAS 
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Solving Second Order ODEs, 

Two Non-analytical Methods Revisited 

Michel BEAUDIN and Gilles PICARD 

Ecole de technologie superieure 
1100, Notre-Dame street West, Montreal (Quebec) Canada 
micheLbeaudin@etsmtl.ca 
gilles. picard@etsmtLca 


ABSTRACT 


Mathematics can still be taught without using a CAS and this is probably the case 
in most schools and universities. Although CAS and technology are often used by 
instructors to demonstrate or illustrate mathematical concepts, they are rarely used 
by students. When we consider our mathematics curriculum, “Differential Equa- 
tions” is one course that we firmly believe can and should benefit from the use of 
CAS. In this talk, we will report how our ODE course has evolved, as our engi- 
neering students have access to technology (Voyage 200 symbolic calculator) in 
the classroom at all times. This talk will show examples of what students still do 
by hand and what CAS allows us to do now to enrich the learning experience. 

We will consider the series solutions of a second order equation with variable coef- 
ficients and numerical solution of the first order equivalent system. Many text- 
books do not show the relation between these two subjects. With a CAS on every 
desk, we can ask students to compare results obtained with both methods. Of 
course, technology is a must to support this. As teachers, we still want our students 
to be able to do some specific computations manually. For example, they have to 
find the recurrence formula by hand for the coefficients of the series solution. 
However, we also want students to be able to compute, with some accuracy, the 
value of the solution at a certain point using partial sums or even graph this ap- 
proximate solution. Then, converting the same equation into a first order system, 
they can plot the numerical generated curve obtained by the built-in RK method in 
the Voyage 200 or create a table of values for the approximate solution. 
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Here are summarized the results to some of the former questions for both problems: 


eql (a= 1) 

eql (a= 2) 

a) c 0 = 4, q = 1, 

(n 2 + 2/i-4) 

c , = - — —c (n > 0). 

n+2 4(/i + l)(n + 2) 

a) c 0 = 8, q = -2, 

(in 2 -3n)c 4c 2 
c , = V ’ (n> 1). 

9n(n + 1) 

c) y(l) = 6.88 (first 5 non vanishing terms). 

c) y( 2) = 0.14 (first 5 non vanishing terms). 

when using more terms, the answer 

when using more terms, the answer 

seems to stabilize to 6.89. 

seems to stabilize to 0.84. 


Figure 2 Results for our two ODEs 


Figure 2 indicates that if we only use the 5 first non vanishing terms, the value ofy(l) will be 
a good estimate in the case of eq 1 while the value of y( 2) will be very bad in the case of eq2. With 
technology, this can be easily corrected. Because of the poor estimate (0.14) obtained in eql , we will 
show how Voyage 200 can help for this equation. 

First, the Y Editor of the Voyage 200, in sequence graphing mode, requires a sequence de- 
fined in the form of “u\(n) =” (see figure 3 a). This means that the students have to rewrite the recur- 
rence formula using the change n —> n - 1 , thus obtaining 

( 2 (« - 1) 2 - 3 (» - l )) c „_, - 4 c ,_ 3 ( 2 « : - 7 » + 5 ) c n _] - 4 c ,_ 3 

9{n - 1) n 9 n(n- 1) 

In Voyage 200, initial values should be given in reverse order (see figure 3 a) and, in this par- 
ticular case where the recurrence formula is of order 3, we need to compute by hand the value of c 2 . 

Recalling that c n = 0 when n < 0, we find 

= (2-2 2 -7-2 + 5)(-2)-4-(0) = 2 = 1 
Cz_ 9-2-1 ~ 18 _ 9 


(a) 

Figure 3 Voyage 200 Y Editor in “sequence” graphic mode (a) and (b) showing a table of values. 
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The coefficients of the series solution are now available through this special function u\(n ) , 
as can be seen in figure 3 b. Partial sums for finding a stable value fory(2) are now computed: 



Figure 4 Function definition and evaluation 


Figure 4 shows us that an appropriate 2 decimal value fory(2) should be 0.84. Using only the 
first 5 non vanishing terms (that is s(4) because c 0 ,c 19 ...,c 4 are all different from 0), we would con- 
clude that y(2) = 0.137174... 

Another way of generating a partial sum for the series solution, without using the graphic en- 
vironment, is to use the “when( , , )” command which is equivalent to a “if then else “ instruction. 
Let’s create a conditional function a(n ) giving us, for n = 0, 1, 2, 3--- the desired coefficients. For the 
above example, the function must give us a value of 0 if n < 0 , 8 if n = 0 , -2 if n = 1 and 
(2 n 2 -ln + 5\a n _ x -4a w _ 3 

if n = 2, 3, 4, 5- • • . Figure 5 a below shows such a function created using 

9n ( n - 1) 

four cascading when commands. 



(a) (b) 

Figure 5 A function to generate the coefficients, and the partial sum of the series solution. 


Looking at figure 5 b, we see one advantage of this approach. The coefficients are in exact 
mode instead of floating point values as in figure 3 b. On the other hand this approach is much slower 
if we need to calculate more precision since evaluating a(n) requires recursion back to 0 for each 
coefficient instead of relying on previous value already calculated as with our previous u\(n) function. 
But the user will get the same values as those obtained using the “sequence” graphic mode (see figure 
6 ). 
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Figure 6 Estimating y( 2) with the a(n) function. 


The Voyage 200 possesses also a robust adaptive RK method [2], so students can do the same 
problem a second time, using another method! This is often not seen in a Differential Equations cour- 
se. The Runge-Kutta methods are usually shown for first order ODE’s. To solve a second order equa- 
tion we need to show students how to transform this equation in a set of first order equations. Of cour- 
se, once this is done, the CAS calculator will do all the computing though students have to learn how 
to control the environment. 

First, let’s look at converting a second order ODE into a system of first order ODE. Every se- 
cond order equation which can be re-written in the form y" = F(x,y,y') can be transformed in a ca- 
nonical way in a system of first order equations. In figure 7 b, we see that the independent variable 
must be t instead of x. Furthermore, the dependant variables must be y 1, y2, y3,. . . In fact, this graphic 
mode is designed to solve system of first order differential equations. 



(a) (b) 

Figure 7 The “differential equations” graphing mode 


Recalling that eql is eql : (2x - 9 )y” -y - 4 xy = 0, y(0) = 8, y'(0) = -2 , we set x = t , then 
y = yl and y f = yY = y2 , we can convert the second order ODE eql into a first order ODE system 


systeql: 


syst _eq 2 


yv = yi 


y 2 ' 


4 ty\ + y2 
2t — 9 


jl(0) = 8 
y2(0) = -2 


We enter these values in the Voyage 200 DE graphing mode (see figure 8 a). One can notice 
the particular way the initial conditions are entered ( £0 = 0 and yi\ = 8 for yl(0) = 8 ). In figure 5 we 

see that we have chosen a graph window with x or £-axis going from -1 to 5 and the y-axis going from 
-5 to 10. 
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F7W 


|Y FiT Y Fi YFH YFEtY FbT Y = " Y~ 

■ |Zqqn |Edit_| y- |fl 1 1 |sty 1 e h .. | 


JLDTl- 

to=o. 

pyl 1 =y2 
yi 1=8 

2 ,_ 4iyl +y2 


yi2= -2 
y3'=ff 
yi3= 
y4 1 = 

■ ■ ' a — 


2 ■ t - 9 


Iv3 1 <t >= 


Bh 


F iH[' fl UTD 


^Tz^r,] 1 

to=o. 
tnax=5. 
tstep=. 1 
tplot=0. 
xn i n= " 1 . 
xnax=5. 
xscl=l . 
ynin= -5. 
ynax=10. 
yscl=l . 
ncur-ues=C 
diftol’FW 
f 1 rlr-p^.=VT 


mm 

2 FiliP liUTD PE 


(a) 


(b) 


Figure 8 Preparing to estimate y( 2) using RK method for syst_eq2 


The step for this example will be 0.1 though this being an adaptive step implementation of a 
Runge-Kutta method the step could be smaller to ensure that the global error will not exceed a value of 
diftol = 0.001 when going from tO = 0 to tmax = 5. In the FI menu, one must also set the “Fields” op- 
tion to “FLDOFF” to prevent the calculator from trying to plot a slope field, option we use when 
studying with them graphical solution of first order equations. Figure 9 shows the curve computed 
with this method and a table of values generated. 








yl 






8. 






.5 

8.9558 






1. 

5.5568 






1.5 

3.5894 






2. 

. 83787 






2.5 

-1.811 






3. 

-3.026 






3.5 

-1.448 






t=0. 

mm 

Rflp flUTD 

QE 1 


(b) 


Figure 9 Solution and table of values obtained using RK method for syst_eq2 


The value 0.83787 was obtained using a step size of 0.1 with a tolerance of 0.001. Changing 
the tolerance to 0.0001 would provide y(2) = 0.83991 . We can be confident with the value 0.84. Note 
also that the graph “stops” at x = 4.5 even though we asked to plot a solution up to t = 5! This method 
being adaptive and wanting to ensure a global tolerance, the calculator isn’t capable of calculating a 
sufficiently small step, thus stopping calculations at the singularity of this equation which is 9/2. Cool! 


4. Conclusion 

This approach — using a mix of paper and pencil and technology — has many advantages. 
One among them is the fact that students are able to “discover” their CAS calculator when they have 
to do this kind of problem using 2 different graphic windows. We have to admit that the RK method 
is only used to check the results obtained by the power series method: we don’t study RK in details 
(only Euler) but this can be done in a separated numerical analysis course. Another advantage of us- 
ing the sequence graphing mode is contained in the following remark: students attending this ODE 
course had followed, earlier, a calculus course where they studied infinite series and convergence. 
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This is a nice topic where we can talk again of these results. They know that if a series of real num- 

oo 

bers y] a n converge, then we must have lima^ = 0. The table of values in figure 3b) shows this but, 

n = 0 

in order to evaluate a series solution of an ODE at some point, “how long will it take” if we want to 
be confident with a partial sum? Experimenting with the CAS Voyage 200 gave us a satisfying an- 
swer. 

Some authors, especially in Applied Differential Equations, will mention that efficient nu- 
merical approximation methods and access to computer and more powerful calculators could question 
the need for power series method in a first course in differential equations, in particular for solutions 
about a regular point. Some of them will even omit this subject [3]. We still think it’s a good idea to 
link all these mathematic topics together and of course, a partial sum of a series solution (a polyno- 
mial) is still a great object to manipulate algebraically. 
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Also, in the writing of this article, we took a look at some exams or exercises given by our 
colleagues to their students. Special thanks to our colleague Chantal Trottier. 


When I came home from Malaga and a wonderful travel through Andalusia I 
studied Michel’s and Gilles’ paper. I had my own ideas for demonstrating the 
convergence behaviour of the power series approximation: 


I set up a table in the Voyage 200 DATA-MATRIX Editor: 



■ 2 a 


2 


1 MAIN 

^^3 hi 1 1 1 

Mi 


CTpiot F W,up' 

feel llHeaderlf 

Tl T FbT Y~7? -I 

:alc|Uiil|stat| 

DATA 

n 


I(c2<n>> 


cl 

lc2 

c3 

14 

15 

16 
17 
IS 
19 
28 

13.00GGGG 

msmsaiH 


14. GQGQGQ 

. QGQG14 

. S4G1G7 

15. GQGGGG 

. QGQG26 

. S4G133 

16. GGGGGG 

. GGGG13 

. S4G146 

17. GGGGGG 

. GGGGG5 

. S4G15G 

IS. GGGGGG 

. GGGGG2 

. S4G152 

19. GQGQGQ 


. S4G152 




Ifir20c3= -840] 

L52401055GG 1 

i iii ijj -itw^^gaE 


This was in 2010. Now we have 2012 and I’d like to transfer the procedure to TI-NspireCAS 
first and then to DERIVE, too. Josef 
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I defined the recurrence equation in the Graphs Application. Take care that the initial 
elements are given in the right order! 





? ■ 



n 

u1(n): T T 





C(2*n ft 2-/*n+S) , u.. 




0. ... 

8.00000000 * 




1 

-2.00000000 




2.... 

0.11111111 




2.... 

*0.58847737 




A.... 

0-02303429 




5 ... 

0.00031243 




$.... 

0. 00875869 




7.... 

0.00038633 

0.2 

. 

* 


0.00014621 

7. M 


1 


y.„. 

-0.00003026 




10-j 

-0 00000992 


m 


11 

-0-00000230 


uiU ( 2 " 3 

7 iW-s] ul(n-l) 4- ul in j) 


12-| 

■0-00000030 





13 

■0-00000003 




14. . 

6. 51 3536C4e-1 0 


, , lta-n a -7- n+r>] ■ ul in- 1 )- a ■ ul fn-s) 



IS... 

, 7.96916C64e-10 


u 1 1 nj 

9 n- fn-l) 



16 

. 1.SB1S4/fc7c-10 



1 



17 

. 3.617332MB-11 


Initial 7inra:«9.*2. 






9 

— 

1S..H 

. 5.7616B1771-12 jU 


nstep*l 



8 

FT 


The spreadsheet is the appropriate application for showing the convergence behavior. The 
value for a can be changed in cell A1. The partial sums are given in column E. 



B alpha 

■n 

c coefficients 

® summands E approx_values 

■ 

♦ 


=seq(k.k,0.40) 



=cumulativesum(summands) 


1 

2 

0 

8.00000000 

8.00000000 

8.00000000 


2 


1 

-2.00000000 

-4.0000000... 

4.00000000 


3 


2 

0.11111111 

0.44444444 

4.44444444 


4 


3 

-0.58847737 

-4.70781 89 J 

-0.26337449 


5 


4 

0.02503429 

0.40054870 

0.13717421 


6 


5 

0.00031245 

0.00999848 

0.14717269 


7 


6 

0.00875869 

0.56055592 

0.70772861 


8 


7 

0.00098633 

0.12624995 

0.83397856 


9 


8 

0.00014821 

0.03794155 

0.87192011 


10 


9 

-0.00003028 

-0.0155030.] 

0.85641711 


11 


10 

-0.00000992 

-0.01 01 553 J 

0.84626178 


12 


11 

-0.00000230 

-0.0047140... 

0.841 54771 


13 


12 

-0.00000030 

-0.001 241 0.J 

0.84030664 


14 


13 

-0.00000003 

-0.00021 40.J 

0.84009259 


15 


14 

8.51355604E-... 

0.00001395 

0.84010653 


16 


15 

7.9891 8064E-... 

0.00002618 

0.84013271 


17 


16 

1 .981 84762e _ ... 

0.00001299 

0.84014570 



The numerical solution of the DE is represented in another Graphs page: 





m 

$.07 ' 



x Utl ,yl_0l(x);= T 





iullerpDlate(!<.de1 . x_0l . 



\ 


1 ,20000000 

4. 6242451 6 - 





1 .30000000 

4.41479142 





1 .40000000 

3.9761 9710 





1 .50000000 

3. S091 2514 





1 .6000MCC 

3.01494933 





1 .70000000 

2.49651 316 





1 .60000000 

1.957 10006 





1 .90000000 

1.40221995 





2. 00000000 

0.83767473 


I 



2.10000000 

0.27124442 




2.20000000 

'0. 2879 796/ 

1 


0.2 

- 4 

2-30000000 

■0,82955150 




\ 

2-4000OOOO 

1.34161 598 




■ 

2-50000000 

1.81114244 





2.60000000 

'2.22415058 



YU 

(2-00000000,0.83767473), X- 2.00000000 

2.70000000 

-?. 50008409 



2-00000000 

*2.82261 870 





1 Ju- 

□ y2 = 



2-00000000 

2.97965072 


S-jr-9 


3.00000000 

3.02624504 - 



o -2 )je£ 

*d 

1 0. «7674727 13432 

l«T4 
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In DERIVE I wanted to do a little bit more. I intended to do all the necessary calculations 
CAS supported: 

n 

#1: y(x) := c -x 

n 

#2: (2-x - 9)*y''(x) - y'(x) - 4-x-yCx) = 0 

n - 2 3 

#3: - x -C4-x + n-x- (3 - 2-n) + 9-n-Cn - l))-c = 0 

n 


n + 1 n - 1 

#4: 4-x -c + n-x -(3 - 2-n)-c 

n n 


n - 2 

9-n-x -(n - l)-c = 0 

n 


I omitted the summation symbols for finding the recurrence equation for the coefficients. 

I performed substitions in the summands of #4 in order to have a common power x n in all of 
them as follows: n — > n - 1 in the 1 st expression, n — ► n + 1 in the 2 nd one and n — ► n + 2 in 
the third one. The rewritten sum is given in #8 where I can factor out x n \ 

n 

#5: 4-x -c 

n - 1 
n 

#6: x -(n + 1)-C1 - 2-n)-c 

n + 1 
n 

#7: 9-x -Cn + l)-(n + 2)-c 

n + 2 

n n n 

#8: 4-x -c + x -Cn + 1)-(1 - 2-n)-c + 9-x -Cn + l)-(n + 2)-c 

n-1 n + 1 n + 2 


#10 gives the recurrence equation - but not in its most comfortable form. 

n 

#9: x -(9*(n + l)-(n + 2)-c + 4-c + Cn + 1)-(1 - 2-n)-c ) 

n + 2 n-1 n + 1 

9-(n + l)-(n + 2)-c + 4-c + Cn + 1)-(1 - 2-n)-c 

#10: n+2 n-1 n+1 


The greatest index should be n. 

So I substitute n — ► n - 2 and obtain #12. 


4-c + Cl - n)-(2-n - 5)-c + 9-n-Cn - l)-c 

#12: n - 3 n-1 n 


Solving this equation for c n needs a little trick. 


S0LVE(4-c + Cl - n)-C2-n - 5)-c + 9-n-Cn - l)-x_, x_) 

#13: n - 3 n-1 

4-c + (1 - n)-(2-n - 5)-c 

n - 3 n-1 

#14: x_ = 

9-n-Cl - n) 

4-c + Cl - n)-C2-n - 5)-c 

n - 3 n-1 

#15: c = 

n 9-n-Cl - n) 


It is not necessary now to calculate c 2 manually! I define the function for the coeffi- 
cients, cf(n), recursively. 
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#27: TAYL0R_0DE2 


v + 4-x-t 


x, t, v, 0, 8, -2, 5 


^ 2-x - 9 


5 4 3 2 

41 -x 73- x 143- x x 

#28: + + 2-x + 8 

131220 2916 243 9 


4828 


#29: 


32805 

#30: 0.1471726870 


#29 is the result of substituting 2 in #28, and #30 is its numerical approximation. 


#31: TAYL0R_0DE2 


v + 4-x-t 
2-x - 9 
10 


x, t, v, 0, 8, -2, 10 


86064547 -x 


#32: - 


9 8 7 6 5 4 

1751803 -x 529297- x 6989 -x 62063 -x 41 -x 73 -x 

+ + + + + 


8678218953600 57854793024 3571283520 7085880 7085880 131220 2916 


#33: Tim TAYL0R_0DE2 

x-»2 


v + 4-x-t 


2-x - 9 


-, x, t, v, 0, 8, -2, 10 


#34: 


57375352024 
67798585575 
#35: 0.8462617846 


#32 is truncated. 

Finally I apply DERIVE s implemented RK-method: 


#36: [yl :=, y2 :r] 


#37: 


RK 


y2, 


4-t-yl + y2 


2-t - 9 


, [t, yl, y2] , [0, 8, -2], 0.005, 400 u[l, 2] 



Tracing the graph of the integral curve gives the value y(2) « 0.84015. 
The respective part of the function table is presented below. 


r / 

4-t-yl + y2 


RK 

y2, 

, [t, yl, y2] , [0, 3, -2], 0.005, 405 


! 2-t - 9 J 

JJ 


[397 403] 


1.93 

0.9535140960 

1.935 

0.9251352125 

1.99 

0 . S96S479S33 

1.995 

0.3635034306 

2 

10 . 34015253591 

2.005 

0.3117964394 


2.01 


0.7334361907 
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Points in a Circle 

Roland Schroder 


Grid points in a circle 


Given is a circle with its center in (0,0) and its 
radius r. We will count all points with integer 
coordinates lying within the circle or on its 
circumference. 

What is the relation between the radius of the 
circle and the number of these points? 


We produce a matrix which contains a “1” for each point which is counted and a “0” for all 
points which are lying outside of the circle. 



F(r ) := 

If a A 2 + b A 2 < r A 2 
# 1 : 1 
0 

#2: G(r) := VECTOR (VECT0R(F(r) , a , - r , r ) , b, -r, r) 

#3: G(4) = 


Nice pattern!! 


0 0 0 
0 0 1 
Oil 
Oil 
111 
Oil 
Oil 
0 0 1 
0 0 0 


0 10 
111 
111 
111 
111 
111 
111 
111 
0 10 


0 0 0 
10 0 
110 
110 
111 
110 
110 
10 0 
0 0 0 


We add all the Ones from G(r) in order to obtain the number P(r) of the grid points within the 
circle with radius r. 


#4: P(r) := X(KG(r))) 

#5: P(4) = 43 


Next step is to create a value table for the 
relation radius r^ number of points P(r): 

#6: VECT0R([r, P(r)], r, 20) 


What is the equation of the graph 
corresponding with this relation? 



200 


10 


20 
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Possible ways to answer might be: 

(1 ) One could assume that it is parabola of form f(x) = a x 2 and then try to find the 
apprpriate a by installing a slider. 



Grid points in the unit circle 

We keep the radius constant r= 1 and increase the number of points by decreasing the 
distance between the points, i.e. we refine the grid: 

#13: H(x) VECTO R (VECTOR (F(l), a , -1, 1, x), b, -1, 1. x) 

For x = 0.1 we get: 



Counting works now as follows: 

#15: count(x) := X(2(H(x))) 

#16: countCO.l) = 317 
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Each point counted is assigned a small square having the point as its center and the side 
length x. The sum of all areas of all these squares is given by: 

2 

#17: area(x) := x -countO) 

#13: area(O.l) =3.17 

#19: areaCO.Ol) = 3.1417 


The value table for function area(x) can easily be created, either using the VECTOR 
command or applying the TABLE function: 


f 


' 0.1 v 

> 

/ 

0.1 

r 0.1 v 

> 

VECTOR 

n, area. 


, n, 10 

TABLE 

, area 


, n, 1, 15 

s 

I 

^ n .j . 

j 

V 

n 

w. n j , 

J 


1 3.17 

2 3.1425 

3 3.134444444 

4 3.140625 

5 3.138 

6 3.135833333 

7 3.137346933 

8 3.13765625 

9 3.141353024 

_ 10 3.1417 


1 

0.1 

3.17 

2 

0.05 

3.1425 

3 

0.03333333333 

3 . 134444444 

4 

0.025 

3.140625 

5 

0.02 

3.133 

6 

0 . 01666666666 

3.135333333 

7 

0.01423571423 

3.137346933 

3 

0.0125 

3.13765625 

9 

0 . 01111111111 

3.141353024 

10 

0.01 

3.1417 

11 

0.009090909090 

3.133925619 

12 

0.003333333333 

3.140625 

13 

0.007692307692 

3.140650337 

14 

0.00714235714 2 

3.139234693 

. 15 

0.006666666666 

3.141377777 


Additional comments of the editor: 

I find this a very nice investigation. It could be followed by the best known Monte Carlo 
method for approxmaing n. 

I wanted to reproduce Roland’s investigation with TI-NspireCAS and was not able to do it in 
the same way because I am missing one of my favourite DERIVE commands: VECTOR. 

So I developped two short functions p(r) and ct(x) and further area(x). 

The Lists & Spreadsheet application and the Graphs Application (including a slider bar) were 
used. Instead of DERIVE’s FIT function (which forced the form of the approximating parabola 
as a-x 2 1 applied a statistics tool, the quadratic regression. 

See the respective screen shots on the next page, Josef 
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The Calculator App together with the two functions: 


p(s) 

81 

seq (£>(&), A, 0,10) 

{ 1,5,13,29,49,81,113,149,197,253,317 } 

4>i) 

317 

o 

d 

31417 

area[x)\=x^ ■ c/(x) 

Done 

area{ O.l) 

3.170000 

area[ O.Ol) 

3.141700 


Define p(r)= 
Func 

Local s,kj 
s:= 0 

For k-r,r 
For j,-r,r 


EndFor 
EndFor 
EndFunc 


a 


s: =s+when(/r 2 +/ 2 <r 2 1 , o) 


"ct" stored successfully 


Define ct(x)= 
Func 

Local s,kj 
s := 0 

For k,-l,l,x 
For 


s: =s+when(& 2 +/ 2 <l,l,oJ 


EndFor 

EndFor 

EndFunc 


5/7 


Part of the spreadsheet together with the scatter diagram and the approximating parabola 
followed by the output of the quadratic regression and the list of the approximated areas of 
the unit circle. 



A n 

B pts 

c distanceBareas E 

l l 

=seq(kj 




=QuadReg( l n,'pts,1 ): 

0 

1 


Title 

Quadratic Regressi... 

1 

5 

0.100000 

3.1 70000 RegEqn 

a*x A 2+b*x+c 

2 

13 

0.050000 

3.142500 a 

3.147892 

3 

29 

0.033333 

3.132222 b 

-0.517393 

4 

49 

0.025000 

3.140625 c 

1.815249 

5 

81 

0.020000 

3.138000 R 2 

0.999979 

6 

113 

0.016667 

3.133611 Resid 

{-0.81524926691,0.... 

7 

149 

0.014286 

3.136939 



p 36 


Nelson Urrego: Shor’s Quantum Algorithm 


D-N-L#86 


USING DERIVE TO SIMULATE THE BASIC STEPS OF 
SHOR’S QUANTUM ALGORITHM 
Nelson Urrego 

Head of the Department of Mathematics 
Gimnasio La Montana 
Bogota, Colombia 

nelsonurrego@glm.edu.co 


Abstract 

It is known that, at the moment, there is no efficient classic algorithm to factor large 
integers. However, in 1994 Peter Shor presented a procedure, which exploits 
quantum parallelism, and is able to efficiently factor large integers. The present work 
describes the basic steps of factoring Shor’s algorithm, and how it can be simulated, 
in a simple way (for small integers), with a CAS like DERIVE. 


Introduction 

The problem of factoring large integers is generally considered as intractable from the Classic 
Theory of Computation. The goal is to find a number k that divides certain number n exactly. 
If k= 1 or k = n, we say that A: is a trivial factor of n. 

As it is known, several cryptosystems (like the famous RSA) base their security on the difficult 
computing problem for factoring large integers. At the present, the best classic algorithms 
to factor integers can execute the task over integers up to 150 digits with a quite demanding 
effort. The running time of these algorithms increases super polynomialy as the integer to be 
factored grows. One of the best algorithms for factoring (The Number Field Sieve) requires 

exp(@(n 1/3 log 273 n)j operations, to factor an /7-bit integer. (Nielsen) 

In 1994 Peter Shor discovered a quantum algorithm that shows that a quantum computer 
can solve the factoring problem exponentially faster than the best known algorithms of its 
counter part classic. With a quantum computer, the task can be accomplished by using 
0(n 2 log n log log «) operations, which represents an important theoretical result in Computer 
Science. (Shor) 


In the first part of this work, we discuss a particular function, from Modular Arithmetic, that 
allows an approach to solve the factoring problem. 
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In the second part we briefly describe the most important stages of Shor’s algorithm and how 
some of them can be simulated with DERIVE. It is clear that some important tools like 
quantum parallelism can not be simulated in a classic computer, but we can visualize the most 
important parts of the procedure with an easy example to factor n = 15 and see how the 
algorithm works. 


A trick from Number Theory 

Consider the function f n (a ) = x mod«, where n is the number to factor and x is a random 
integer co-prime with n. The function f„(a) is periodic, which means that as you increase 
the argument a = 0, 1,2, ..., k, their values /„( 0), /„( 1), /„( 2), ..., f n (k) turn in some pattern 

for a given x. The number of values between such pattern is called the period or the order of 
x modulo n, and is denoted by r. 

This implies that: 


If r is even we have: 


x r = lmodw 


(x W2 ) 2 = lmodra 
(x W2 -l)(x W2 + l) = 0mod« 


( 1 ) 


( 2 ) 


and (x W2 - l)(x W2 +1) is a factor of n. 


Now unless that x W2 = ±lmod«, at least one of the terms on the left hand side in (2) has a 
non- trivial common factor with n. Hence by computing gcd(x' 2 - 1 ,n) and gcd(x r/2 + 1 ,n) we 
have a good chance to find a factor of n. The problem is that, at the present, there is no 
efficient algorithm to find r on a classic computer. 


Steps of Shor’s Algorithm 

Shor’s algorithm is an example of a quantum algorithm, this is, a procedure which can be 
carried out by using a quantum system. 


The main goal of Shor’s algorithm is to obtain efficiently the period r of the function 
f n (a ) = x" mod n. In order to accomplish this, the algorithm exploits quantum parallelism. Let 
n be the number to be factored, for instance n= 15. 
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Step 1 

Pick a number q such that In 2 <q< 'in 2 . This guarantees that q is not so large. 

This can be programmed with DERIVE: 

2 2 
q := 2-n + RANDOM (3 -n ) 

q := 243 
Step 2 

Pick a number x which is co-prime with n. 

x := L00P(x := RANDOM (n) , IF(GCD(x, n) = 1, RETURN x, x := x + 1)) 
x :z 13 

Step 3 

Create a single quantum memory register which we divide into two sets called reg 1 and regl. 
The combined register is represented by | regl,reg2). 

Step 4 

Load reg\ with integers between 0 and q - 1 and regl with zeros 
regl := VECTORS + 1, a, -1, q - 1) 
reg2 := VECT0R(0-m, m, 0, q - 1) 

The combined state is represented by | y/) = — ^=^|g,0). 


Step 5 

Apply * a mod n to each number in regl and write the result in regl. The combined register 


becomes \a,x° mod«). 

\Jq a=0 


n := 15 


reg2 := VECTOR ( 
reg2 := [1, 13, 
7, 1, 13, 4, 
1, 13, 4, 7, 
13, 4, 7, 1, 
4, 7, 1, 13, 


(M0D(x , n) , ; 
, 4, 7, 1, 13 
, 7, 1, 13, 4 
, 1, 13, 4, 7 
, 13, 4, 7, 1 
, 4, 7, 1, 13 


a, 0, q - 1) 

, 4, 7, 1, 13 
, 7, 1, 13, 4 
, 1, 13, 4, 7 
, 13, 4, 7, 1 
, 4, 7, 1, 13 


4, 7, 1, 13 
7, 1, 13, 4 
1, 13, 4, 7 
13, 4, 7, 1 
4, 7, 1, 13 


4, 

7, 

1, 

13 

7, 

1, 

13, 

, 4 

1, 

13 

, 

, 7 

13 

, 4 

, 7, 

, 1 

4, 

7, 

1, 

13 


7, 1, 13, 4 


7, 1, 13, 4 


7, 1, 13, 4 


7, 1, 13, 4 


7, 1, 13, 4] 
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In this sample, for n = 15, q = 243 and x = 13, it is easy to conclude that the Period r = 4, 
However we are interested in describing the complete algorithm in a general form and 
estimate the period by using the discrete Fourier transform. 


Step 6 

Measure the state in register reg2 and obtain certain k. In quantum computation this has the 
consequence of projecting the state of regl into a superposition of all values a' such that 
x a mod n = k. 

RANDOM (q - 1) 

MeasureReg2 := M0D(x , n) 


The complete state of the quantum register is |^) = where 


A = \a':x a modn = k\ and A represents the number of elements in A. 


projectRegl := VECTOR (I F(MOD(x , n) = MeasureReg2, b, 0), b, 0, q - 1) 


ketft) := 


1 q - 1 

reglKet := ■ I ket(t) 

fo t=0 

b 

z := DIM (SELECT (MGD(> , n) = MeasureReg2, b, regl)) 

b 

reglmeasure := SELECT (M0D(x , n) = MeasureRegZ, b, regl) 


1 z 

■ I ket(ELEMENT(reglmeasure , d)) 

v'z d=l 


Step 7 

Compute the discrete Fourier transform of the projected state in regl : 


Inia'ct q 


V 1 

rZ-rS- 

| cigA yj c = 0 


f 2 Tuia'c! q 


FourierRegl := VECTOR 


■ ELEMENT(reglmeasure , d)-c/q 


VECT0R([s - 1, | ELEMENT (FourierRegl , s) | ] , s, 1, q - 1) 
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0 

0.4885520725 ' 

1 

0.02309639933 

2 

0.0228938444 

3 

0.02255814073 

239 

0.02209210406 

240 

0.02255814073 

241 

0.0228938444 

1 

0.02309639933 . 


1.2568 


Hi . 0054 


0.75408 


0.50272 


0.25136 


q/r 


-H- 


q/r 




q/r 




q/r 


S ^ 


37.755 75.51 113.27 151.02 


188.78 


226.53 264.29 


-0.25136 


0.50272 


Regl [tf : jc fl mod/; = 13] 


Figure 1 is the graph of register 1 which contains the set {a: x a mod n = 13} in the horizontal 
axis and ln(pr(measure a)) in the vertical axis. 

Here the discrete Fourier transform projects the state of regl into a superposition 



1 1 q ~ l 


e 


Ima'c! q 


c, 


k). 


Step 8 

By measuring the state of regl we obtain samples of the discrete Fourier transform. The result 
is certain number t which is multiple of q/r (Figure 1). We repeat 0(ln q) times the steps 1 - 7 
to estimate period r. 

In this case r = 4, therefore by computing gcd(13 4 2 - 1,15) and gcd(13 4/2 + 1,15) we obtain the 
desired factors 3 and 5 respectively. 
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Conclusion 

Recently Shor’s algorithm has been implemeted by using new technologies like the employ of 
spins in liquid solution. There is still a long way to cover before we can see a functional 
quantum computer, but several theoretical results have risen in the last decade, involving the 
discovery of quantum algorithms, which allow us approaching problems, which are not 
feasible for classic computers. Such results, along with other important outcomes in quantum 
cryptography (Urrego) make the quantum information processing a fascinating field to 
explore and to study. 
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I found some interesting websites, Josef. 

Related websites: 

http://arxiv.org/pdf/quant-ph/9508027v2.pdf (Peter Shor’s original paper) 

http://de.wikipedia.org/wiki/Shor-Algorithmus (in German) 

http://www.worldcat.org/title/quantum-computation-and-quantum-information/oclc/43641333 

(Nielsen & Chang’s book) 

http://en.wikipedia.org/wiki/Shor%27s algorithm 
http://katzgraber.org/teaching/FS08/files/herrigel.pdf 

http://www.cs.washington.edu/education/courses/cse599d/06wi/lecturenotesl 1 .pdf 

http://www.newscientist.com/blog/technology/20Q7/09/how-quantum-computer-factorises- 
numbers.html (A nice short explanation) 


About Quantum Computers 

http://en.wikipedia.org/wiki/Quantum computer 

http://computer.howstuffworks.com/quantum-computerl.htm (How QuantumCompters Work) 
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Samples from Dietmar Oertel’s exciting Taylor Series Gallery 
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In the last DNL we mentioned the TRIBONACCI-Sequence (Cubus Simus contribution). 
I invited writing a function/program to produce this sequence. 

With DERIVE we can do this performing a RECURSION or performing an ITERATION : 


tribo(n) :z 
If n z 0 
0 

If n z 1 
1 

If n = 2 
1 

If n = 3 
2 

tribofn - 1) + tribo(n - 2) + tribo(n - 3) 
VECTOR (tribo(k), k, 1, 10) 

[1, 1, 2, 4, 7, 13, 24, 44, 81, 149] 


a „ = «„-! + a n-2 + «„-3 M » ^ 4 

a A = a 2 = 1, c/ 3 = 2 


Recursion (above) and iteration (below) 


# 1 : 

# 2 : 

#3: 

#4: 

#5: 


tnb_seq(n) := (ITERATES(|~w , 


w , 
3 


w + w + w , 
1 2 3_ 


w, 


[1, 1, 2], n - l))ul 


tnb_seq(15) 

[1, 1, 2, 4, 7, IB, 24, 44, SI, 149, 274, 504, 927, 1705, 3136] 
tnb_el(n) := (trib_seq(n)) 

n 


trib_el (15) = 3136 


With TI-Nspire we could write a program / function or we use the appropriate mode in the Graphs 
Application and produce the function table. 


6.74 ; 
1 • 

k y 

m 

m 

■ m m 

X 

n 

u1(n):= T 

T 

ul (n-1 )+u1 (n-2)+u1 (n-3) 


1 

1 

s 

2. 

37 

1 


2 


4. 

4 


lll(tf)=Ul(/7-l)+U: 

l(/7— 2)+Ul(tf— 3) 

5. 

7 


6. 

13 




77 

24 


Ul(77)=Ul(n-l)+Ul(/7-2)+Ul(/7-3)| 
| Initial Terms- =1,1,2 
l<n<99 nstep= 1 


8. 

44 


9 

81 

H 

1 


MM 


\ wj(io) 


0 

149. f 

seq(i/7(4*,l,20) {l.,l.,2., 4*7,13,24,44*8 1,1 49, 274, 504*927, 1705,3 136,5768, 10609, 19513,35890,66012} 


The sequence can be called in the Calculator App using the seq -command. 
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Dear Josef, 

I cannot find an example of a function with an arbitrary number of arguments. 
Thank you very much. 

Best regards. 

Robert Setif 


I don’t know if you are looking for an easy example or a more sophisticated one. 
I add an example using DERIVE’s FOO-function. 


# 1 : 


# 2 : 


#3: 


foo v := [AVERAGE (v) , VARIANCE(v)] 


i 

o 

£ 

8) = 

1 1 
I- 1 

UJ UJ 

1 

37 ' 
3 _ 


i 

o 

£ 

S, 10 

II 

l n 
1 

' 3S 293 

_ 3 

1 

i- 1 
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Dear Josef, 

Thank you very much for your fast answer. 

If you have time, please would you send me several or at least one example of a more sophisticated 
one. 

Best regards 

It would be great if anybody could send another example! 


Dear Josef, 

I need to construct a matrix of all zeros (32 x 32) and then go 
through the point plot that I make and depending on the location 
of the point, update the matrix by adding 1 to the corresponding 
position in the matrix. Nothing seems to happen. The matrix 
stays all zeros. I tried the small example to try and debug my 
larger program. I attached the small file. I tried several strategies 
with the same non-result. What am I missing? Is there a utility 
file that I need to load? 


See my answer: 



' 1 

2 

3 

4 ' 


A := 

_ 3 

6 

7 

9 _ 


A 

:= A 


+ 

1 


2,3 


2,3 





' 1 

2 

3 

4 ' 


A = 

_ 3 

6 

7 

9 _ 


REPLACE (A 

■+ 

■ 1 , 

A , 3) 



2, 

3 


2 


' 1 

2 

3 

4 ' 


A = 

_ 3 

6 

7 

9 _ 



#1: A := 


12 3 4 

3 6 7 9 


A := A + 1 
#2: 2,3 2,3 

Now simplify 


Or more directly: 

Enter expression #2 again and press the —button 

in the entry line 

for immediately simplifying: 


A := s 
#3: 2,3 


#3: 


A := 9 
2,3 


The element has changed as expected: 


A = 


1 2 3 4 1 


# 6 : 


A = 


" 1 2 3 4 " 

_ 3 6 § 9 _ 


#4: 


_ 3 6 B 9 _ 


